Variant profiling of colorectal adenomas from three patients of two families with MSH3-related adenomatous polyposis.

The spectrum of somatic genetic variation in colorectal adenomas caused by biallelic pathogenic germline variants in the MSH3 gene, was comprehensively analysed to characterise mutational signatures and identify potential driver genes and pathways of MSH3-related tumourigenesis. Three patients from two families with MSH3-associated polyposis were included. Whole exome sequencing of nine adenomas and matched normal tissue was performed. The amount of somatic variants in the MSH3-deficient adenomas and the pattern of single nucleotide variants (SNVs) was similar to sporadic adenomas, whereas the fraction of small insertions/deletions (indels) (21-42% of all small variants) was significantly higher. Interestingly, pathogenic somatic APC variants were found in all but one adenoma. The vast majority (12/13) of these were di-, tetra-, or penta-base pair (bp) deletions. The fraction of APC indels was significantly higher than that reported in patients with familial adenomatous polyposis (FAP) (p < 0.01) or in sporadic adenomas (p < 0.0001). In MSH3-deficient adenomas, the occurrence of APC indels in a repetitive sequence context was significantly higher than in FAP patients (p < 0.01). In addition, the MSH3-deficient adenomas harboured one to five (recurrent) somatic variants in 13 established or candidate driver genes for early colorectal carcinogenesis, including ACVR2A and ARID genes. Our data suggest that MSH3-related colorectal carcinogenesis seems to follow the classical APC-driven pathway. In line with the specific function of MSH3 in the mismatch repair (MMR) system, we identified a characteristic APC mutational pattern in MSH3-deficient adenomas, and confirmed further driver genes for colorectal tumourigenesis.

Genotype Pattern Mining for Pairs of Interacting Variants Underlying Digenic Traits.

Some genetic diseases ("digenic traits") are due to the interaction between two DNA variants, which presumably reflects biochemical interactions. For example, certain forms of Retinitis Pigmentosa, a type of blindness, occur in the presence of two mutant variants, one each in the ROM1 and RDS genes, while the occurrence of only one such variant results in a normal phenotype. Detecting variant pairs underlying digenic traits by standard genetic methods is difficult and is downright impossible when individual variants alone have minimal effects. Frequent pattern mining (FPM) methods are known to detect patterns of items. We make use of FPM approaches to find pairs of genotypes (from different variants) that can discriminate between cases and controls. Our method is based on genotype patterns of length two, and permutation testing allows assigning p-values to genotype patterns, where the null hypothesis refers to equal pattern frequencies in cases and controls. We compare different interaction search approaches and their properties on the basis of published datasets. Our implementation of FPM to case-control studies is freely available.

A distinct APC pathogenic germline variant identified in a southern Thai family with familial adenomatous polyposis.

BACKGROUND: Familial adenomatous polyposis (FAP) is caused by pathogenic germline variants in the APC gene. To date, multiple pathogenic variants in coding regions, splice sites, and deep intronic regions have been revealed. However, there are still pathogenic variants that remain unidentified. METHODS: Twenty-nine primer pairs flanking exons 2-16 (i.e., coding exons 1-15) of APC and their exon-intron junctions were used for germline pathogenic variant screening in Southern Thai patients with familial adenomatous polyposis (FAP). Transcription analysis was performed to confirm the pathogenicity of a splice site deletion of intron 10. Family members were interviewed for clinical histories. Blood samples were collected from 18 family members for a segregation study. Subsequently, clinical data of affected members were collected from the hospital databases. RESULTS: We found a distinct heterozygous 16-bp deletion at the splice donor site of intron 10 leading to a skipping of exon 10 which was confirmed by transcript analysis (APC: c 1312 + 4_1312 + 19del, r.934_1312del). Predictive testing for the pathogenic APC variant in 18 of the proband's family members (ten healthy and eight affected) from three generations showed the same heterozygous germline pathogenic variant in eight affected adult members (15-62 years old) and two children (7 and 10 years old). Seven of the ten carriers of the disease-causing variant had undergone colonoscopy, and colonic polyps were found in all cases, which confirmed the segregation of the inherited pathogenic variant. The phenotypic spectrum was found to vary within the family; and some affected family members exhibited extracolonic manifestations. CONCLUSIONS: To our knowledge, the pathogenic APC variant, c.1312 + 4_1312 + 19del, r.934_1312del, has not previously been reported. This study is one of the few reports describing the phenotypic consequences of a pathogenic APC variant in a high number of affected family members.

Shared genomic segment analysis with equivalence testing.

An important aspect of disease gene mapping is replication, that is, a putative finding in one group of individuals is confirmed in another set of individuals. As it can happen by chance that individuals share an estimated disease position, we developed a statistical approach to determine the p-value for multiple individuals or families to share a possibly small number of candidate susceptibility variants. Here, we focus on candidate variants for dominant traits that have been obtained by our previously developed heterozygosity analysis, and we are testing the sharing of candidate variants obtained for different individuals. Our approach allows for multiple pathogenic variants in a gene to contribute to disease, and for estimated disease variant positions to be imprecise. Statistically, the method developed here falls into the category of equivalence testing, where the classical null and alternative hypotheses of homogeneity and heterogeneity are reversed. The null hypothesis situation is created by permuting genomic locations of variants for one individual after another. We applied our methodology to the ALSPAC data set of 1,927 whole-genome sequenced individuals, where some individuals carry a pathogenic variant for the BRCA1 gene, but no two individuals carry the same variant. Our shared genomic segment analysis found significant evidence for BRCA1 pathogenic variants within ±5 kb of a given DNA variant.

Heterozygosity mapping for human dominant trait variants.

Homozygosity mapping is a well-known technique to identify runs of homozygous variants that are likely to harbor genes responsible for autosomal recessive disease, but a comparable method for autosomal dominant traits has been lacking. We developed an approach to map dominant disease genes based on heterozygosity frequencies of sequence variants in the immediate vicinity of a dominant trait. We demonstrate through theoretical analysis that DNA variants surrounding an inherited dominant disease variant tend to have increased heterozygosity compared with variants elsewhere in the genome. We confirm existence of this phenomenon in sequence data with known dominant pathogenic variants obtained on family members and in unrelated population controls. A computer-based approach to estimating empirical significance levels associated with our test statistics shows genome-wide p-values smaller than 0.05 for many but not all of the individuals carrying a pathogenic variant.

Copy number variation analysis and targeted NGS in 77 families with suspected Lynch syndrome reveals novel potential causative genes.

In many families with suspected Lynch syndrome (LS), no germline mutation in the causative mismatch repair (MMR) genes is detected during routine diagnostics. To identify novel causative genes for LS, the present study investigated 77 unrelated, mutation-negative patients with clinically suspected LS and a loss of MSH2 in tumor tissue. An analysis for genomic copy number variants (CNV) was performed, with subsequent next generation sequencing (NGS) of selected candidate genes in a subgroup of the cohort. Genomic DNA was genotyped using Illumina's HumanOmniExpress Bead Array. After quality control and filtering, 25 deletions and 16 duplications encompassing 73 genes were identified in 28 patients. No recurrent CNV was detected, and none of the CNVs affected the regulatory regions of MSH2. A total of 49 candidate genes from genomic regions implicated by the present CNV analysis and 30 known or assumed risk genes for colorectal cancer (CRC) were then sequenced in a subset of 38 patients using a customized NGS gene panel and Sanger sequencing. Single nucleotide variants were identified in 14 candidate genes from the CNV analysis. The most promising of these candidate genes were: (i) PRKCA, PRKDC, and MCM4, as a functional relation to MSH2 is predicted by network analysis, and (ii) CSMD1, as this is commonly mutated in CRC. Furthermore, six patients harbored POLE variants outside the exonuclease domain, suggesting that these might be implicated in hereditary CRC. Analyses in larger cohorts of suspected LS patients recruited via international collaborations are warranted to verify the present findings.

Exome sequencing characterizes the somatic mutation spectrum of early serrated lesions in a patient with serrated polyposis syndrome (SPS).

BACKGROUND: Serrated or Hyperplastic Polyposis Syndrome (SPS, HPS) is a yet poorly defined colorectal cancer (CRC) predisposition characterised by the occurrence of multiple and/or large serrated polyps throughout the colon. A serrated polyp-CRC sequence (serrated pathway) of CRC formation has been postulated, however, to date only few molecular signatures of serrated neoplasia (BRAF, KRAS, RNF43 mutations, CpG Island Methylation, MSI) have been described in a subset of SPS patients and neither the etiology of the syndrome nor the distinct genetic alterations during tumorigenesis have been identified. METHODS: To identify somatic point mutations in potential novel candidate genes of SPS-associated lesions and the involved pathways we performed exome sequencing of eleven early serrated polyps obtained from a 41 year-old female patient with clinically confirmed SPS. For data filtering and analysis, standard pipelines were used. Somatic mutations were identified by comparison with leukocyte DNA and were validated by Sanger sequencing. RESULTS: The BRAF p.V600E or KRAS p.G12D mutation was identified in six polyps (~50%) and not found in polyps from the distal colon. In addition, we found seven unique rare somatic alterations of seven different genes in four serrated tumours, all of which are missense variants. The variant in ABI3BP and CATSPERB are predicted to be deleterious. No established cancer gene or candidate genes related to serrated tumorigenesis were affected. CONCLUSIONS: Somatic mutations seem to be rare events in early hyperplastic and serrated lesions of SPS patients. Neither frequently affected genes nor enrichment of specific pathways were observed. Thus, other alterations such as non-coding variants or epigenetic changes might be the major driving force of tumour progression in SPS.

Exome Sequencing Identifies Biallelic MSH3 Germline Mutations as a Recessive Subtype of Colorectal Adenomatous Polyposis.

In ∼30% of families affected by colorectal adenomatous polyposis, no germline mutations have been identified in the previously implicated genes APC, MUTYH, POLE, POLD1, and NTHL1, although a hereditary etiology is likely. To uncover further genes with high-penetrance causative mutations, we performed exome sequencing of leukocyte DNA from 102 unrelated individuals with unexplained adenomatous polyposis. We identified two unrelated individuals with differing compound-heterozygous loss-of-function (LoF) germline mutations in the mismatch-repair gene MSH3. The impact of the MSH3 mutations (c.1148delA, c.2319-1G>A, c.2760delC, and c.3001-2A>C) was indicated at the RNA and protein levels. Analysis of the diseased individuals' tumor tissue demonstrated high microsatellite instability of di- and tetranucleotides (EMAST), and immunohistochemical staining illustrated a complete loss of nuclear MSH3 in normal and tumor tissue, confirming the LoF effect and causal relevance of the mutations. The pedigrees, genotypes, and frequency of MSH3 mutations in the general population are consistent with an autosomal-recessive mode of inheritance. Both index persons have an affected sibling carrying the same mutations. The tumor spectrum in these four persons comprised colorectal and duodenal adenomas, colorectal cancer, gastric cancer, and an early-onset astrocytoma. Additionally, we detected one unrelated individual with biallelic PMS2 germline mutations, representing constitutional mismatch-repair deficiency. Potentially causative variants in 14 more candidate genes identified in 26 other individuals require further workup. In the present study, we identified biallelic germline MSH3 mutations in individuals with a suspected hereditary tumor syndrome. Our data suggest that MSH3 mutations represent an additional recessive subtype of colorectal adenomatous polyposis.

Exome sequencing identifies potential novel candidate genes in patients with unexplained colorectal adenomatous polyposis.

In up to 30% of patients with colorectal adenomatous polyposis, no germline mutation in the known genes APC, causing familial adenomatous polyposis, MUTYH, causing MUTYH-associated polyposis, and POLE or POLD1, causing Polymerase-Proofreading-associated polyposis can be identified, although a hereditary etiology is likely. To uncover new causative genes, exome sequencing was performed using DNA from leukocytes and a total of 12 colorectal adenomas from seven unrelated patients with unexplained sporadic adenomatous polyposis. For data analysis and variant filtering, an established bioinformatics pipeline including in-house tools was applied. Variants were filtered for rare truncating point mutations and copy-number variants assuming a dominant, recessive, or tumor suppressor model of inheritance. Subsequently, targeted sequence analysis of the most promising candidate genes was performed in a validation cohort of 191 unrelated patients. All relevant variants were validated by Sanger sequencing. The analysis of exome sequencing data resulted in the identification of rare loss-of-function germline mutations in three promising candidate genes (DSC2, PIEZO1, ZSWIM7). In the validation cohort, further variants predicted to be pathogenic were identified in DSC2 and PIEZO1. According to the somatic mutation spectra, the adenomas in this patient cohort follow the classical pathways of colorectal tumorigenesis. The present study identified three candidate genes which might represent rare causes for a predisposition to colorectal adenoma formation. Especially PIEZO1 (FAM38A) and ZSWIM7 (SWS1) warrant further exploration. To evaluate the clinical relevance of these genes, investigation of larger patient cohorts and functional studies are required.

Low-level APC mutational mosaicism is the underlying cause in a substantial fraction of unexplained colorectal adenomatous polyposis cases.

BACKGROUND: In 30-50% of patients with colorectal adenomatous polyposis, no germline mutation in the known genes APC, causing familial adenomatous polyposis, MUTYH, causing MUTYH-associated polyposis, or POLE or POLD1, causing polymerase-proofreading-associated polyposis can be identified, although a hereditary aetiology is likely. This study aimed to explore the impact of APC mutational mosaicism in unexplained polyposis. METHODS: To comprehensively screen for somatic low-level APC mosaicism, high-coverage next-generation sequencing of the APC gene was performed using DNA from leucocytes and a total of 53 colorectal tumours from 20 unrelated patients with unexplained sporadic adenomatous polyposis. APC mosaicism was assumed if the same loss-of-function APC mutation was present in ≥ 2 anatomically separated colorectal adenomas/carcinomas per patient. All mutations were validated using diverse methods. RESULTS: In 25% (5/20) of patients, somatic mosaicism of a pathogenic APC mutation was identified as underlying cause of the disease. In 2/5 cases, the mosaic level in leucocyte DNA was slightly below the sensitivity threshold of Sanger sequencing; while in 3/5 cases, the allelic fraction was either very low (0.1-1%) or no mutations were detectable. The majority of mosaic mutations were located outside the somatic mutation cluster region of the gene. CONCLUSIONS: The present data indicate a high prevalence of pathogenic mosaic APC mutations below the detection thresholds of routine diagnostics in adenomatous polyposis, even if high-coverage sequencing of leucocyte DNA alone is taken into account. This has important implications for both routine work-up and strategies to identify new causative genes in this patient group.

Leveling the Playing Field in Homozygosity Mapping Using Map Distances.

Studies of linkage disequilibrium (LD) and its variation in the genome are of central importance for understanding evolutionary history, population structure, and selective sweeps. Extreme forms of the latter may result in runs of homozygosity (ROH). In human gene mapping, long ROHs are the basis for homozygosity mapping (HM) with length measured in terms of Mb (106 base pairs physical distance). LD varies greatly over the human genome so that long ROHs tend to occur preferentially in regions of high LD and ROHs of the same length in different regions are not strictly comparable. Thus, in human gene mapping, LD appears as a confounder that needs to be taken into account in the interpretation of ROHs. The effect of varying LD can be mitigated by working on a scale of centimorgans (cM, genetic distance) instead of Mb. We demonstrate this effect for HapMap 3 data on chromosome 19 and show examples with different ROH lengths depending on whether physical or genetic lengths are used. These results suggest that HM should preferably be done on genetic rather than physical distances.

Genome-wide CNV analysis in 221 unrelated patients and targeted high-throughput sequencing reveal novel causative candidate genes for colorectal adenomatous polyposis.

To uncover novel causative genes in patients with unexplained adenomatous polyposis, a model disease for colorectal cancer, we performed a genome-wide analysis of germline copy number variants (CNV) in a large, well characterized APC and MUTYH mutation negative patient cohort followed by a targeted next generation sequencing (NGS) approach. Genomic DNA from 221 unrelated German patients was genotyped on high-resolution SNP arrays. Putative CNVs were filtered according to stringent criteria, compared with those of 531 population-based German controls, and validated by qPCR. Candidate genes were prioritized using in silico, expression, and segregation analyses, data mining and enrichment analyses of genes and pathways. In 27% of the 221 unrelated patients, a total of 77 protein coding genes displayed rare, nonrecurrent, germline CNVs. The set included 26 candidates with molecular and cellular functions related to tumorigenesis. Targeted high-throughput sequencing found truncating point mutations in 12% (10/77) of the prioritized genes. No clear evidence was found for autosomal recessive subtypes. Six patients had potentially causative mutations in more than one of the 26 genes. Combined with data from recent studies of early-onset colorectal and breast cancer, recurrent potential loss-of-function alterations were detected in CNTN6, FOCAD (KIAA1797), HSPH1, KIF26B, MCM3AP, YBEY and in three genes from the ARHGAP family. In the canonical Wnt pathway oncogene CTNNB1 (ß-catenin), two potential gain-of-function mutations were found. In conclusion, the present study identified a group of rarely affected genes which are likely to predispose to colorectal adenoma formation and confirmed previously published candidates for tumor predisposition as etiologically relevant.

Frequency and phenotypic spectrum of germline mutations in POLE and seven other polymerase genes in 266 patients with colorectal adenomas and carcinomas.

In a number of families with colorectal adenomatous polyposis or suspected Lynch syndrome/HNPCC, no germline alteration in the APC, MUTYH, or mismatch repair (MMR) genes are found. Missense mutations in the polymerase genes POLE and POLD1 have recently been identified as rare cause of multiple colorectal adenomas and carcinomas, a condition termed polymerase proofreading-associated polyposis (PPAP). The aim of the present study was to evaluate the clinical relevance and phenotypic spectrum of polymerase germline mutations. Therefore, targeted sequencing of the polymerase genes POLD1, POLD2, POLD3, POLD4, POLE, POLE2, POLE3 and POLE4 was performed in 266 unrelated patients with polyposis or fulfilled Amsterdam criteria. The POLE mutation c.1270C>G;p.Leu424Val was detected in four unrelated patients. The mutation was present in 1.5% (4/266) of all patients, 4% (3/77) of all familial cases and 7% (2/30) of familial polyposis cases. The colorectal phenotype in 14 affected individuals ranged from typical adenomatous polyposis to a HNPCC phenotype, with high intrafamilial variability. Multiple colorectal carcinomas and duodenal adenomas were common, and one case of duodenal carcinoma was reported. Additionally, various extraintestinal lesions were evident. Nine further putative pathogenic variants were identified. The most promising was c.1306C>T;p.Pro436Ser in POLE. In conclusion, a PPAP was identified in a substantial number of polyposis and familial colorectal cancer patients. Screening for polymerase proofreading mutations should therefore be considered, particularly in unexplained familial cases. The present study broadens the phenotypic spectrum of PPAP to duodenal adenomas and carcinomas, and identified novel, potentially pathogenic variants in four polymerase genes.

MUTYH-associated polyposis (MAP): evidence for the origin of the common European mutations p.Tyr179Cys and p.Gly396Asp by founder events.

MUTYH-associated polyposis (MAP) is an autosomal recessive adenomatous polyposis caused by biallelic germline mutations of the base-excision-repair gene MUTYH. In MAP patients of European origin, the combined allele frequency of the mutations p.Tyr179Cys and p.Gly396Asp ranges between 50 and 82%, while these mutations have not been identified in Far Eastern Asian populations, supporting the hypothesis that a founder effect has occurred at some point in European history. To investigate the natural history of the two common European MUTYH alleles, we genotyped six gene-flanking microsatellite markers in 80 unrelated Italian and German MAP patients segregating one or both mutations and calculated their age in generations (g) by using DMLE+2.2 software. Three distinct common haplotypes, one for p.Tyr179Cys and two for p.Gly396Asp, were identified. Estimated mutation ages were 305 g (95% CS: 271-418) for p.Tyr179Cys and 350 g (95% CS: 313-435) for p.Gly396Asp. These results provide evidence for strong founder effects and suggest that the p.Tyr179Cys and p.Gly396Asp mutations derive from ancestors who lived between 5-8 thousand years and 6-9 thousand years B.C., respectively.

Deep intronic APC mutations explain a substantial proportion of patients with familial or early-onset adenomatous polyposis.

To uncover pathogenic deep intronic variants in patients with colorectal adenomatous polyposis, in whom no germline mutation in the APC or MUTYH genes can be identified by routine diagnostics, we performed a systematic APC messenger RNA analysis in 125 unrelated mutation-negative cases. Overall, we identified aberrant transcripts in 8% of the patients (familial cases 30%; early-onset manifestation 21%). In eight of them, two different out-of-frame pseudoexons were found consisting of a 167-bp insertion from intron 4 in five families with a shared founder haplotype and a 83-bp insertion from intron 10 in three patients. The pseudoexon formation was caused by three different heterozygous germline mutations, which are supposed to activate cryptic splice sites. In conclusion, a few deep intronic mutations contribute substantially to the APC mutation spectrum. Complementary DNA analysis and/or target sequencing of intronic regions should be considered as an additional mutation discovery approach in polyposis patients.